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equation plus an under-determined ODE (that is, a semi-explicit DAE system of dif- 
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where, for any integer i s.t. 1 < i < n, /j is a polynomial in the variables X := {xi, . . . , x n } 
and in their jth (1 < j < e^) time derivatives A^- 7 ') := {xi^', . . . ,x n ^}, with coefficients 
in a differential field IK of characteristic 0. 

One of the main invariants of DAE systems is their differentiation index. There 
are several definitions of differentiation indices not all completely equivalent (see for in- 
stance [H El HH EH [291 [371 EH SU [121 S31 SH S9] ) , but in every case it represents a measure 
of the implicitness of the given system in a fixed coordinates set. For instance, for first 
order equations, differentiation indices provide bounds for the number of total derivatives 
of the system needed in order to obtain in the same set of coordinates an explicit ode 
system which is verified by all the solutions of the original system (see (H Definition 2.2.2]). 

Since explicitness is strongly related to the existence of classical solutions, a differ- 
entiation index should also bound the number of derivatives needed in order to obtain 
existence and uniqueness theorems (see |40 1 141 1 l4"3]). From the point of view of numerical 
resolution, it is desirable for the DAE to have an index as small as possible. As shown in [H 
§2.5.3], for first order systems a reduction of the index can be achieved by differentiating 
the algebraic constraints, but the numerical solution of the resulting system do not satisfy 
necessarily the original equations. 

Main contributions. In this article we address the index reduction problem for an 
ubiquitous class of quasi- regular DAE systems (see Section I2.3() . We show that any of 
these systems is generically equivalent to a related (in a non intrinsic way) first order 
DAE system S with a particular structure. This new system consists of a single purely 
algebraic (polynomial) equation plus an under-determined ode (see Definition [3l page fill) . 
Indeed, £ is a semi- explicit DAE system in the usual sense (see for instance [H Section 
1.2]) with differentiation index 1 (see Proposition 1101 page I15j) . It is a well-known fact 
that this class of systems can be handled successfully by means of numerical methods 
(see [381 EH E]). 

The index reduction problem has already been considered in several previous articles 
(see for instance [151 [T6| [26l [34"] ) . The techniques applied in these works are based on the 
computation of sufficiently many successive derivatives of the original equations combined 
with rewriting procedures relying on the Implicit Function Theorem, elimination of critical 
equations, introduction of dummy derivatives, etc. 

Our approach also makes use of the computation of successive derivatives, as many as 
the index, but, unlike the methods mentioned above, we deal with the system of all these 
new equations in a purely algebraic way. This new system defines an algebraic variety in a 
suitable jet space and we parametrize this variety by means of the points of a hypersurface. 
This construction, originally introduced by Kronecker, is known as a geometric resolution 
(see [El [471 [13] an d references therein). In order to keep track of the differential structure, 
we use the parametrizations to construct a vector field over the hypersurface defining the 
semi-explicit DAE system E. A result of the same flavor (i.e. an univariate differential 
equation plus parameterizations of the variables) may be given by means of the notion of 
primitive element of an extension of differential fields. This construction, due to J. Ritt 
([46], see also [IE]), is known as a resolvent representation of the system £ (see [6] \7\ [TT] 
for effective versions of it). 

With respect to the known index reduction approaches, our method is symbolic and, 
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in some sense, automatic: it does not make use of the Implicit Function Theorem as it 
is the case in [15] or [16] and it does not rely on any smart choice of ad-hoc equations 
as in [26] . Moreover, the construction can be done algorithmically within an admissible 
complexity by applying well-known techniques from computer algebra (see [47} 130]). 

The number of variables of our semi-explicit system is the order of the differential ideal 
associated to the input system plus one and it is always lower than those involved in the 
index reduction methods of previous papers. In the first order case, where it is easy to 
compare, this number is at most n + 1 and in the general case, it is bounded by the Jacobi 
number of the system (see [20j [Ml EH] and [10] ) . 

A further advantage of our method is that it preserves the constraints of the initial 
conditions of the original system and then we do not need to compute constants of inte- 
gration. 

Outline. The paper is organized as follows: in Section [2] the basic notions needed 
throughout the article are introduced. The core of the paper is in Section [3] where the 
semi-explicit system E is constructed. In Section U] we study the relation between the so- 
lutions of both systems. Finally, two appendices are included: the first one contains some 
Bertini type results from commutative algebra we need and the second one is devoted to 
existence and uniqueness theorems for DAE systems. 

2 Preliminaries 

In this section, we introduce some notations used throughout this paper and we recall 
some basic definitions from elementary (differential) algebraic geometry for the reader 
convenience. Furthermore, we discuss the assumptions on the systems considered and 
some results concerning the differentiation index and the order of these systems. 

2.1 Basic notions and notations 

Let Kbe a characteristic zero field equipped with a derivation 5. For instance K = Q, R 
or C with 5 = 0, or K = Q(i) with the usual derivation 5(t) = 1, etc. 

As in the Introduction, for any set of X := {xi, . . . ,x n } of n (differential) indetermi- 
nates over K, we denote by x%j> the j-th successive formal derivative of the variable x% 
(following Newton's notation, its first derivative is also denoted by x{) and we use the 
following notations: 




and 




The derivation 5 can be extended to a derivation in the polynomial ring j £ No] 

as follows: for any differential polynomial q in KjX^), j 6 No] the following classical 
recursive relations hold for the successive total derivatives of q: 





for j >1 
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where 6\qV 1 ^| denotes the polynomial obtained from qU ^ by applying the derivative 5 
to all its coefficients. The (non-Noetherian) polynomial ring K[X^, j € No] with this 
derivation is denoted by K{a;i, . . . , x n } (or simply K{JT}) and is called the ring of differ- 
ential polynomials. 

Given a finite set of (differential) polynomials Q := {q±, . . . , q u } in we write [Q] 

to denote the smallest ideal of stable under differentiation, i.e. the smallest ideal 

containing qi, . . . , q v and all their derivatives of arbitrary order. The ideal [Q] is called the 
differential ideal generated by Q. Furthermore, for every integer j, we extend our previous 
notations as follows : 

Q® := {qi®,...,q„®} and Q® := [q, Q, , . . . , Q^). 

Let us introduce also some notions concerning elementary algebraic geometry. 

Let Y := {y\, . . . , y m } be (algebraic) indeterminates over the field K; we write K[V] to 
denote the polynomial ring in m variables over K. Let K be a fixed algebraic closure of IK. 
Given some polynomials pi, . . . ,p p in K[Y"], the set {y £ IK , pi(y) = • • • = p p (y) = 0} 
is called an algebraic variety definable over K (or simply a variety if IK is clear from the 
context). The affine space K is endowed with a topology (the so-called Zariski topology) 
where the closed sets are exactly the algebraic varieties definable over K. We denote this 
topological space by A" 1 . The space A m is a Noetherian space and then every closed set 
is an irredundant union of a finite number of irreducible closed sets. 

Given a subvariety V of A m we denote by /(V) the ideal in K[Y] of all the polynomials 
that vanish on V and by K[V] := K[y]/7(V) the coordinate ring of V. 



2.2 The considered system — Primality assumption 



In this section we recall some notations concerning the DAE systems considered in this 
paper. Then, we explicit a natural primality assumption necessary in the sequel. 

Let n denote a fixed non-negative integer. Throughout the paper we consider DAE 
systems of the following type: 



E := 



h(x,x,...,x^ 



fjx,x,...,x^ 



0, 



0, 



where, for every 1 < i < n, fi is a polynomial in the variables X and the derivatives X^\ 
with 1 < j < ei] the coefficients of these polynomials are in the field K. Each non-negative 
integer denotes the maximal derivation order appearing in the polynomial /j. We 
write e := max{ej} for the maximal derivation order that occurs in £ and we assume 
that e is greater or equal to 1. As done previously, we use the following notations: 

F:= {h,. ..,/„}, F' J : {./V ./V'} and F® := {f, F, . . . , F®}. 

Let [F] C be the differential ideal generated by the polynomials F. We introduce 

also the following auxiliary (Noetherian) polynomial rings and ideals: for every j in No, R"^ 
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denotes the polynomial ring K[XW] and pr^F the ideal in R0'- 1+e ) generated by the total 
derivatives of the defining equations up to order j — 1, namely pr^F := (i 7 ^ -1 !) (this ideal 
is usually known as the (j — l)th prolongation ideal). We set pr^F := (0) by definition. 

For i = 0, . . . , n in N and for every integer j, we will assume that the ideals generated 
by the polynomials i 7 ^" 1 !, f\^\ . . . , ffi' are all prime ideals in their respective rings. In 
particular, the differential ideal [F] is a prime differential ideal in the ring K{A}. 

2.3 Quasi-regular DAE systems and prime complete intersection 

In this section, we establish a relationship between the notion of quasi-regularity of a 
differential system and an algebraic property — complete intersection — that is required by 
the geometric elimination algorithm used in this paper. The notion of quasi-regularity 
appears implicitly in [22] in order to generalize a conjecture of Janet to non-linear systems. 

Definition 1 — Let T be a DAE system given in the ring IK{A} by differential poly- 
nomials Q := {q±, . . . ,q u } of order bounded by a nonegative integer e. Let p be a prime 
differential ideal containing [Q]. We say that T is quasi-regular at p if for every integer j 
in No, the Jacobian matrix of the polynomials with respect to the set of variables X^ +e ^ 
has full row rank over the domain R(J+ e )/(R(-?' +e ) n p). We say that T is quasi-regular if it 
is quasi-regular at any minimal prime differential ideal containing [Q\. j 

For the systems £ considered in this paper, since the ideal [F] is assumed to be prime, 
the quasi-regularity of £ is equivalent to say that for each integer j, the Jacobian ma- 
trix of the polynomials F^ with respect to the set of variables X^" 1 " 6 ! has full row rank 
over the domain R(J+ e )/(R(-J+ e ) n [F]). This condition can be easily rephrased in terms 
of Kahler differentials (as in Johnson's original work [22]) saying that the set of differen- 
tials {dfi ij \ 1 < i < n, j e N } are a K{ X}/[.F]-linearly independent set in the module of 
differentials ^k{x}/k ^k{x} Geometrically, it means that the algebraic variety 

defined by the ideal generated by the (j + l)n polynomials F^ in the ( j + e + l)n-variate 
polynomial ring R( J+e ) is smooth at almost every point of the closed subvariety defined by 
the prime ideal R^+ e ) n [F] . 

Under our assumptions we have the following straightforward consequence: 

Proposition 1 If the system S is quasi-regular, then for i = 0, . . . , n and every j in No, 
the polynomials F^~ l \ f\^\ . . . , fft' form a regular sequence in the ring R^'+ £ ). In par- 
ticular, the prolongation ideals pft'F are prime complete intersection ideals. 

Proof. Let lj, be the ideal of R^ +e ) generated by F^ -1 !, f\^\ . . . , ffi\ From the Jacobian 
Criterion (see [331 §30]) we observe that the quasi-regularity condition implies that the 
ideal L-j is complete intersection and its generators form a regular sequence in the local 
ring (R^' +e ^) , where q denotes R^ +e ^> n [F]. Since the ideal \ji is assumed to be prime 

and it is contained in q, we conclude that the polynomials F^ -1 ^, fi^\ ■ ■ ■ , f%^ form a 
global regular sequence in R^' +e ). ■ 

2.4 Differentiation index — Linearized standpoint 

We introduce here the notion of differentiation index of quasi-regular DAE systems used in 
this paper and establish a relationship between this index and the dimension of Jacobian 
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matrix kernels. We keep the hypotheses on the system £ made on Section [2.21 and, from 
now on, we also assume that S is quasi-regular. 

Consider the following chain C of (prime) ideals in the polynomial ring R^ 6 " 1 ); 

C : = pr^F n R (e ~ 1} C pr^F n R (e_1) C • • • C pr^F n R (e " 1} C • • • C [F]nR< e-1) . 

Since R( e_1 ) is a Noetherian ring, the ideal chain C eventually becomes stationary. Clearly, 
the biggest proper ideal of the chain must be [F] n R^ 6 " 1 ). 

Definition 2 — The differentiation index a of the system £ is the minimum integer j 
at which the chain C becomes stationary; more precisely, 

a := min U G N | pr (j) F n R (e ~ 1} = pr^ j+h ^F n R (e ^ 1} , V/i G n|. 

Clearly we have a = min{j G N | pr^F n R (e ~ 1} = [F] n R^" 1 )}. j 

The differentiation index can also be defined by means of Jacobian matrices related to 
the input system. 

For any positive integers j and h, with h > e — 1, let Zj,h be the Jacobian matrix 
of F( h - e+1 ), . . . ,F( h - e+ J) with respect to the variables X^ h+l \ . . .,X^ h+ ^. Let the inte- 
ger dj t h denote dim(ker( i ^j i /j)) (where t Zj,h denotes the transposed matrix of Zj,h) and 
let do,/i := 0- Under the additional hypothesis that the rank of the matrices Zj,h ls inde- 
pendent of the ring where it is computed (the rank over the rings R^ +h+% ) j '^(i+ h+l - e+1 ) p 
or over Rfr+i+V/dF] n R^+i+V) is the same for any integer i), it can be shown that 
the double sequence djh is in fact independent of h (see [121 Proposition 11]). If we 
write dj := dj hi we have the following alternative characterization of the differentiation 
index (see [101 Theorem 8 and Definition 9]): 

a = min{j G No | dj = dj+i}. 

From this characterization, we deduce the following result (see |1Q(. Theorem 10]): 

Theorem 2 The differentiation index a satisfies: 

a = min{j G N | pr^ +h - e+1 ^F n R (h) = p^ +h ~ e+2 ^F n R (h) } 
= minO G N | p^ +h - e+1 ^F n R (h) = [F] n R (h) } 

for every integer h > e — 1. 

The techniques used in [10] rely on the structure of the Jacobian matrices involved. 
Here we give an alternative proof of the above result for the case h = e based on the 
characteristic set theory (see [23j E5] ) . 

Lemma 1 — If, for some integer j, pr^Fn R( e ) = pr( J+1 )Fn R( e ), then pr^Fn R( e ) = 

[F]nR( £ ). j 

Proof. Let A be an algebraic characteristic set of the prime ideal pr^F n R (e) for some 
orderly ranking on derivatives. From A we extract a minimal chain B as follows: from all 
the polynomials in A with the same leading variable we take the one with the minimal 
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order of derivation in this variable. We claim that B is autoreduced in the differential 
meaning. 

This is equivalent to the fact that, if xf 1 is the leading derivative of some element B 
of B, then this derivative does not appear in some other element. As we use an orderly 
ranking, the e — h first derivatives of B belong to R^ e ^ and, since by assumption pr^ +1 >FH 
R (e) = pr (j) F n R (e) 5 they b e i on g to pr (j) F n R(e)_ Theil) t]ie derivatives xf\ h < £ < e, 

are the leading derivatives of these elements of pr^F n R( e ). These derivatives appear 
with degree 1 and with initial equal to Sb, the separant of B, that does not belong to 
prC?)_p n R( e ). So they are the leading derivatives, with degree 1, of some elements of A, 
and they do not appear in other elements of this characteristic set. Hence our claim. 

So, B is the characteristic set of some prime differential ideal V C [F] (see [2]). Now, 
it is easily seen that all polynomials in pr^F n are reduced to by B, which implies 
that F C V, so that [F] = V, and also that pr^F n R (e) = [F] n R (e) . ■ 

In the last part of this section, we recall in a geometric framework the notion of initial 
conditions associated to a given differential system. 

2.5 Hilbert-Kolchin regularity — Independent variables 

Under our assumptions, the differential dimension of the prime differential ideal [F] is 
(see e.g. |25|). So, following |23} Chapter II, Section 12, Theorem 6], the transcendence 
degree of the fraction field of the domain R^ /(R^) n [F]) over the ground field K becomes 
constant for all i sufficiently big. This constant is a non-negative integer called the order 
of [F] and it is denoted by ord[F]. 

The minimum of the indices jo such that the order of [F] equals the transcendence 
degree of the fraction field of R^/(R W) n [F])) over K for all j > jo is known as the Hilbert- 
Kolchin regularity of the ideal [F]. In our situation, the Hilbert-Kolchin regularity of [F] 
is bounded by e — 1 (see |10l Theorem 12]). 

From the results of the previous subsection, it follows that the differentiation index of 
the system £ is at most en — ord[F] (for more precise bounds, see for instance, [10j). 

Since the fraction fields of the domains R^" 1 ) /(R^" 1 ) n [F]) and R( e ) /(R( e ) n [F]) have 
the same transcendence degree over K, from the canonical inclusion 

R (e-l)/( R (e-l) n ^ ^ R(e)/(RW n [F]), 

we conclude that there exists in X^ e ~^ a subset U of ordfi 7 ] many variables that is a 
transcendence basis of both these fields. Moreover, we may also choose U in such a way 
that Xj G U implies G U for every < £ < h, e.g. U may be chosen as the set 
of derivatives that are not leading derivatives of the algebraic characteristic set A in the 
proof of Lemma [TJ We are going to see in the sequel that this set of variables could be 
considered as initial conditions. 

3 A related vector field over an algebraic hypersurface 

In this section we exhibit a new DAE system S related (in a non intrinsic way) to the 
original one S. This new DAE system has a very particular structure: a single purely 
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algebraic (polynomial) equation q = plus an under-determined ode system (see Defini- 
tion [3] page Q31) • In particular, X is a semi-explicit DAE system in the usual sense (see for 
instance [H Section 1.2]). Moreover, we will prove that the differentiation index of X is 1 
(see Proposition \W\ page fT5|) . 

The polynomial equation g = is obtained by means of a classical, purely algebraic 
procedure known today as the geometric resolution (see Section 13. 2p applied to a suitable 
algebraic variety associated to the input DAE system X and some of its derivatives (see 
Section 13. 1 1) . The differential equations of X are introduced in Section 13.41 

We leave for Section H] the analysis of the relations between the solutions of both DAE 
systems X and X. 

3.1 The prolonged algebraic system and its partial specialization 

We keep the notations and assumptions introduced in Section [2] related to the DAE input 
system X. 

We recall that U denotes a subset of Xl e_1 ] that is a transcendence basis of the fraction 
fields of the domains R^/QF] n R( e_1 )) and R (e) /([^] H R (e) ). Following SectionESsuch 
a basis exists and its cardinality is ord[F]. Recall that a denotes the differentiation index 
of X introduced in Section 12.41 

Proposition 3 1. The variables U as elements of the ring R( e + CT )/p/ CT+1 ).F remain 
algebraically independent over K. 

2. Let W be a subset of X^ +e ^ such that {U, W} is a transcendence basis of the fraction 
field of R( e+ff ) /pM +l ^F. Then every variable in W has order at least e + 1; in other 
words, W is a subset of {X^'; e + l<j<e + a}. 

Proof. Note that Theorem [2] (for i = e) or Lemma Q] imply that the canonical inclusion 
of R (e) in R( e+<T) induces an injective K-algebra morphism R( e ) /{[F] n R (e) ) ^ R (e+<T) /pr^ +1) F. 
In particular, this inclusion preserves K-algebraically free elements and then the state- 
ment (1) follows. In order to prove the second assertion simply observe that U is a transcen- 
dence basis of the fraction field R( e )/([F] n R^ e ^) and then, for every 1 < i < n, {U, x^} is 
an algebraically dependent set modulo [F] D R( e \ and the same holds in R^+^/pr^+^F. 
■ 

Let If be a subset of {X®; e + 1 < j < e + cr} verifying the second assertion in 
Proposition [3] (observe that if a = there are no variables W); since pr^ a+1 ^F is a complete 
intersection prime ideal of the polynomial ring R( e+CT ) (Proposition [1]) , we have that the 
cardinality of {U, W} equals the number of variables of the polynomial ring R( e+<J ) minus 
the number of elements of the regular sequence defining pr ( - cr+1 - ) i ? . In other words: 

card{£7, W} = dim R (e+a) - (a + l)n = (e + a + l)n - (a + l)n = ne. 

Let s be the cardinality of W, that is ne — ord[F]. For any differential polynomial / 
in K{A} and any point W in A s denote by f\y\> the polynomial obtained by replacing in / 
the variables W by the corresponding value W. 
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Proposition 4 There exists a nonempty Zariski open subset of A s such that for any W in 
this set and for all integer i such that 1 < i < a + 1 , the following conditions are satisfied: 

1. The sequence F^lw forO<j<i-l is a reduced regular sequence in 
particular, the ideals pr^>F + (W — W) in R( e+cr ) are radical and complete intersec- 
tion. 

2. No prime component of these ideals p^F + (W — W) contains a nonzero polynomial 
pure in U. 

Proof. The proposition is a consequence of the results given in Appendix [A] The first 
statement follows directly from Theorem [18] and for the second one we apply Corollary [T71 
(remark that the ideals pr^F are supposed to be prime). ■ 

Now we introduce an algebraic variety defined by the prolonged equations of the input 
system T, up to order a followed by a specialization of the variables W. Fix a specialization 
point W in A s and suppose it belongs to the Zariski open set given by Proposition UJ 

Notation 5 Let p/ CT+1 ) F\ w be the ideal spanned by the subset F^\ w ofK[X^ e+ ^ \ W] . 
We denote by V the algebraic (equidimensional) variety in ^ ord [ F \+( 1 + (T ) n defined by the 
ideal pr^- a+l ^ F\yy and by Vi U • • • U Vn the irreducible decomposition ofV. 

Example 1 — The ideal pr( a+1 ^ F\w may actually fail to be prime for all values of W in 

(2) (2)^ 

a dense set, as shown by the following example: x w - x 2 ' =0,x 2 = 0. It is easy to see 
that <T = 2 and we may choose {x\ , x^} as the set W. Then, for an orderly ordering, an 
algebraic characteristic set (in fact, a system of generators) of pr^+^Fjw; is 2x^ — Wi, 

(3) (2) (2) 

x\ , x\ , x 2 , x' 2 , X2- We see that the ideal is prime if and only if Wi/2 is not a square 
in K. Moreover, even in the prime case, the field extension associated to pr( a+1 ^ F\yy is a 
non-trivial algebraic extension of degree 2 of the field associated to the ideal [F] n R^ e \ 
which is 1. On the other hand, we could choose also W as the set {x^\ x^}; in this case 
p r (cr+i)^| w ^ s p r j me anc j [i s associated variety is birational equivalent to VQi 7 ] n R*- 6 )). 

Finding whenever possible, such a choice of W, remains a subject for further investi- 
gations. —Ex H 

Observe that the algebraic variety V is not intrinsically associated to the input DAE 
system because its definition depends on the choice of the transcendence basis U, the 
variables W, and the point W where the variables W are evaluated. 

Let us also remark that the second assertion in Proposition [5] states that the projection 
on the [/-space of any irreducible component Vi is dominant; i.e. the closure of the image 
of Vi by the projection on the variables U is the whole space A 01 "^! or equivalently, the 
natural ring map M,[U] — > K[Vi] is injective. 

The following proposition shows that the identity [F] n R^ e ^ = pr^+^i 71 n R^ e ^ (see The- 
orem [2] and Lemma [T]) remains correct after specialization in a suitable W: 

Proposition 6 Let W in A s chosen as in Proposition ^ Then the identity [F] n R^ e ^ = 
p^ + ^F\ w n R( e ) holds. 
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Proof. Since [F] n R( £ ) = pr^+^F n R( £ ) and the variables W do not appear in the ring R^, 
the ideal [F] n R (e) is included in pr^ +1 ^F\ w n R (e) . On the other hand if p is a primary 
component of pr( a+1 ^ F\w , we have that ~K[U] R( e+CT )/p (because of the choice of W 
verifying Proposition H]). Then p n R^ e ^ is a prime ideal of dimension at least ord[F] con- 
taining [F] n R( e ), which is a prime ideal of dimension ord[-F]. Hence both prime ideals 
are the same. Since the argument holds for any primary component of pr^ +1 ^ F\w> , the 
proposition follows. ■ 

In other words, this proposition says that all differential conditions of order at most e 
induced by the input system can be generated by differentiation of the original equations 
up to order a followed by the specialization W h- > W. 

In particular, if K = R, C, suppose that <p := (ipi, . . . , <p n ) is a classical analytic solution 
of the DAE system £ defined locally in a neighborhood of 0. Then, Proposition [6] implies 
that for any t in R small enough, the complex vector formed by the derivatives up to order e 
of the function <p evaluated at the instant t is a point of the algebraic var 
independently of the choice of the variables U, W and the point W in A s . 

3.2 A parametric geometric resolution of the variety V 

In this Section we introduce the algebraic part of our related semi-explicit DAE system 
£ and certain rational functions that allow us to express solutions of £ from solutions 
of £. To do this, we use a classical tool in effective Algebraic Geometry: the geometric 
resolution of an equidimensional variety. In our case this construction will be applied to 
the algebraic variety V introduced in Notation [5J 

Let us explain informally this well-known notion (for simplicity, we assume IK is alge- 
braically closed and of characteristic 0): suppose that a d-dimensional irreducible affine 
variety V in the m-dimensional ambient space A m is given. Then the field K(V) of ra- 
tional functions over V has transcendence degree d over the ground field K. Therefore, 
there exist d variables, say y\, . . . , yd, such that the extension field . . . , yd) ■— > K(V) 

is finite. These variables are called parametric or free variables. The Primitive Element 
Theorem (see for instance \28\ §V.4, Theorem 4.6]) asserts that there exists an element v 
in K(V) such that K(V) = K(yi, . . . ,yd)[v]; moreover, the element v can be taken as a 
generic K-linear combination of the remaining variables yd+i, ■ ■ ■ , Dm- The minimal poly- 
nomial q of v over K(yi, . . . , yd) defines an irreducible hypersurface T~L := {q = 0} in the 
affine space A d x A 1 (q can be taken with coefficients in the polynomial ring . . . , yd])- 
Since each of the (non-free) variables yd+i, ■ ■ ■ ,y m , as elements of the field K(V), can be 
written as rational functions in the variable v over the field K(yi, . . . , yd), it follows that 
a dense open subset of the irreducible variety V can be rationally parametrized from a 
dense open subset of the hypersurface H. 

The 4-tuple consisting of the parametric set {yi, . . . ,yd}, the element v, its minimal 
polynomial q, and the rational parametrizations is called a parametric geometric resolution 
of the variety V. 

If the variety V is not irreducible but equidimensional, a similar construction can be 
reproduced with suitable changes (see for instance [Ul Section 2]). For instance, there is in 
general no chance that the same choice of the free variables as a subset of the input variables 
works for any component of V. In this case a (generic) linear change of coordinates may 
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be necessary in order to obtain the same parametric set of variables for any irreducible 
component. We point out also that in this case, for any choice of the element v, each 
irreducible component of the induced hypersurface % parametrizes generically one (and 
only one) component of V. In particular, the number of irreducible components of V and % 
is the same. 

In our situation, we consider a parametric geometric resolution for the equidimen- 
sional variety V introduced in Notation From Proposition 0] we observe first that 
the variables U are a parametric set with respect to the equidimensional algebraic vari- 
ety V := Vi U • • • U Vat, since the canonical morphism K[J7] — > K[Vj] is injective and since 
the relations ord[-F] = card U = dim Vi hold for all i = 1, . . . , N. In particular no linear 
change of coordinates is necessary in order to obtain free variables with respect to the 
irreducible components of V. Secondly, the ideal pr( a+1 ^ F\yy is radical and so, it is the 
defining ideal of V. Moreover, it is generated by the regular sequence -F^|>v. 

These facts imply that for each prime ideal /(Vi) associated to defined 
in Notation I(Vi) <S> K(£/) is a O-dimensional prime ideal in the polynomial ring with 
coefficients in ~K(U) and variables A^ 6+<J ] \ {U, W}. Hence the Jacobian determinant of 
the polynomials i^ CT J|w with respect to these variables X^ e+a ^ \ {U, W} does not vanish 
identically over any component Vi. 

Thus the requirements of |47t Section 2.1] are fulfilled and a parametric geometric 

resolution \ U, v, q, | V 1 < i < (1 + a)nj \ exists. Here v is a Q- linear combination of 

the variables X^ e+ ^ \ {U,W}, q the square- free polynomial in K[U, v] of positive degree 
in v defining a hypersurface % in A ord ^] x A 1 , and q' the partial derivative The 
fractions ^ i n K(C7, v) are the parametrizations of the remaining variables. More precisely, 
each pi can be written as 

_ ai(U,v) 

Pl ~ h(u) {2) 

where ai(U,v) in K[{7, v] and bi(U) in ¥L[U] \ {0} are coprime polynomials verifying that 
deg^ ai < deg v q; furthermore, for each variable y in X^- e+a ^ \ {U, W}, there exists an inte- 
ger j such that bj(U)q'y — aj(U,v) vanishes on the variety V. 

We define the total ring of fractions of the variety V in the usual way as the Artinian 
ring K(V) := K(Vi) x • • • x K(Vat) and analogously for T~L. Therefore, from the canonical 
ring inclusions K[C7] M,[H] <—> K[V], by means of the geometric resolution and passing 
to the total ring of fractions, we infer that the relations K(C7) <—t K.(H) = K(V) hold. 
The inverse application K(V) — > K('H) is induced by the parametrization y i— > for y 

in X[ e+<T 1 \{U,W}. 

From a more geometrical point of view, these facts can be stated in the following way. 
Consider the linear map * : A ord [ F H> +1 ) n -»• A ord [ F ] +1 defined as X^ e+ ^ \ {W} ^ (U,v). 
For each irreducible component Vi of V the restriction of ^ to Vi induces an isomor- 
phism between suitable nonempty Zariski open sets of Vi and of the irreducible com- 
ponent %i = ^(Vi) of T~L. In other words, the components Vi and %i are birationally 
equivalent, i.e. the fields of rational functions K(Vj) and K('Hj) are K-isomorphic. 

As shown in |47j . if the polynomials defining V are encoded by straight line programs, 
a parametric geometric resolution can be computed by means of a probabilistic algorithm 
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of bounded complexity in terms of certain parameters. 

In order to estimate the running time of this algorithm in our case, we point out that 
from a straight-line program of length L encoding the input polynomials F, we can obtain 
a straight-line program of length 0(((e + a)en + L)a 2 ^ encoding all the polynomials 
(see |32[ Section 5.2] or jll[ Lemma 21]). Then applying [471 Theorem 2] to our algebraic 
system {i^lyy = 0} and the parametric variables U, we deduce the following complexity 
result concerning the computation of a geometric resolution of V: 

Proposition 7 Assume the input polynomials F are encoded by a straight-line program 
of length L. Then, a parametric geometric resolution of V can be computed through the 
following steps: 

1. take a point U in K ord ^ and compute a geometric resolution of the zeros of the 
system obtained by specializing U = IA in F^|w; 

2. apply a formal Newton lifting process, that requires 

°* (( H^m^lm ) m.<**m.«H 

operations in K, 

where N = ordfF] + (1 + o~)n, L = ((e + a)en + L)a 2 and p stands for the degree of the 
projection V — > A ord [ F l mapping X^- e+a ^ \ W to U ; M(i) denotes the cost of the arith- 
metic operations with univariate polynomials of degrees bounded by i with coefficients in a 
ring — we can take M(i) = 0(i log 2 (i) log log(i)) — and M s (i, j) the cost of j-variate series 
multiplication at precision i — that can be taken less than Oi^M^' 7 )) . 

Remark 1 — It is easy to see that the variables U can be chosen in such a way that for 
every variable Xj £ U, all its previous derivatives xf \ for < I < h, also belong to U . 
In this case the representation given in Proposition [7] allows us to obtain a characteristic 
set of the ideal for some ranking on derivatives, using for instance the method described 
in [9]. j 

We may also consider a particular kind of parametric geometric resolution of the vari- 
ety V that we call a Noether parametric geometric resolution which requires an additional 
condition: the natural morphism K[f7] — > K[V] must be not only injective but also integral 
(i.e. it verifies Noether's Normalization Lemma). This new requirement implies a more 
stable geometric behavior which allows improvements in the algorithms performance. 

We are not able to ensure the existence of a set of variables U for which the variety V is 
in Noether position, but this can be achieved by a (generic) Q-linear change of coordinates 
(see for instance [IT] and [TH]). 

Taking into account that the polynomials -F^yy form a reduced regular sequence, 
we can apply the algorithm presented in [17] to compute a Noether parametric geometric 
resolution of the variety V. This leads to the following complexity result. 
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Proposition 8 Using the previous notations, assume also that the polynomials F have 
degrees bounded by a positive integer d. Then, a Noether parametric geometric resolution 
of V can be computed over the field IK by means of a probabilistic algorithm which runs in 
time 

I / flog 2 (D)l 

O N (NL + N w ) M(Dcf) 2 + M(D) ^ a(2*) 

where D is the maximum of the degrees of the varieties successively defined by the polynomi- 
als F^jw, and, for every positive integer j, a(j) is the cost of the arithmetic operations in 
the quotient R/m 3 , where R is a polynomial ring with coefficients in K in oid[F] variables 
and m is the maximal ideal generated by all the variables (uj denotes the linear algebra 
constant). 




3.3 Computing a, U and W 

Up to now, we have assumed the sets a, U and W to be known a priori. This may often be 
the case for U for obvious physical reasons (for example, one does not need to compute the 
equations of a mechanical system such as the pendulum to know which quantities could 
be arbitrarily chosen). Such an assumption is much harder to justify for W, but we will 
see that suitable sets U and W may be computed with little extra cost. 

According to Theorem [2] or Lemma [1] in order to compute a it is enough to find the 
minimum jo such that pr^) FdR^ = pr^ +1 ) jP(~lR^ e ' (then this minimum jo is cr+1). By the 
primality assumption of these ideals it suffices to compare their dimensions. Following [12\ 
Proposition 2 and Remark 3] or |10t Proposition 6] the dimension of the ideal pr^'FH 
is equal to (e — j + l)n + rank^F^ / 'dX^)i< r< j t e <h<e+j, the rank being computed 
modulo the prime ideal pr^F. Our algorithm computes successively the ranks for the 
ideals pr^'F and it stops when two consecutive dimensions coincide. 

To do this we apply the Kronecker algorithm described in |17l I13j . For every j it 
computes a Noether geometric resolution for the algebraic variety defined by pr^F. This 
geometric resolution allows us to reduce the rank computation modulo pr^'F to a prob- 
abilistic rank computation modulo a principal ideal. Due to the recursive structure of 
Kronecker algorithm, if the equality of the dimensions does not hold, the geometric reso- 
lution already computed can be taken as input for the next step. 

Once the differential index is obtained, then we are able to compute the set of variables 
U and W . This can be done by considering the Jacobian matrix J a of the polynomials 
F^ with respect to the variables X^ e+a \ . . . ,X,X . After a Gauss triangulation of J a , 
the variables indexing columns with no pivot give a transcendence basis modulo pr^^F. 
The set U corresponds to those variables of order at most e — 1 and the set W to the 
remaining ones 

The complexity of this procedure is similar to the one of Proposition [8] for the ad- 
equate parameters, namely: the number of variables and equations; the degrees of the 
intermediate varieties. 

The previous computations might simplify the obtention of a parametric geometric 
resolution of the variety V. This question, and its computational interest, are left to 
further investigations. 
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3.4 An associated vector field over the hypersurface 

In this section we define a vector field on the algebraic hypersurface % defined in A 1+ord [ F ] 
by {q = 0} an d introduced in the previous section. Moreover, we introduce the new first- 
order, quasi-regular system £ having differentiation index 1, whose solutions will enable 
us to obtain solutions of the given system S. 

Consider a parametric geometric resolution of the variety V and let 

(u,v,q,{f,l<i< (l + ^)n}) 

be the parametric variables, the primitive element, its minimal square- free polynomial and 
the parametrizations respectively, as in the previous section. 

The linear map f : A ord [ F l + ( 1+CT ) n -> A 1+ord [ F l defined as \ W -> (U,v) (re- 

call that v is a Q-linear combination of the variables xt e+cr l \ {U, W}) gives, by restric- 
tion, a morphism of algebraic varieties between V and T~L and so, it induces a dual K- 
morphism ^* between the Artinian rings K(%) and K(V). From the properties satisfied 
by the geometric resolution, we have that \&* is an isomorphism of IK-algebras and its 
inverse morphism is defined, by means of the parametrizations, as the dual of the 
(rational, not necessarily polynomial) morphism of algebraic varieties: $ : % — > V defined 
by (U, v) — > (U, ^r|l < i < (1 -+- er)n). Let us observe that both ring morphisms ^f* and 
fix the variables U. 

Since the parametric set U has been chosen as a subset of X^- e ~ 1 ^ , the set U of derivatives 
of U is included in and so, by Proposition [3l the relation U D W = holds. In 
particular, U and U remain invariant after specialization of the variables W at any point W 
in K s . 

Fix a variable vn of the set U (1 < i < ordfi 71 ]). We have: 

(a) If ili is in U, there exists a unique integer h such that 1 < h < ord[F] and Ui = Uh- 

(b) If Ui is not in U, there exists a unique index j such that 1 < j < (1 + a)n and: 

q bj{U) q' 

Definition 3 — Let S be the square DAE system in the ord[-F] + 1 differential un- 
knowns U, v: 

ili — u h = 0, for all ili verifying condition (a) 

bj(U)q' (U,v)ui — dj(U,v) = 0, for all u\ verifying condition (b) 
q(U,v) = 0. 

We denote by F := fi, . . . , A+ rd[F] the polynomials in K[U, U, v] defining the system S 
and by [F] the differential ideal generated by them in BC{f7, v}. j 

Note that S is a semi-explicit DAE system, i.e. an explicit under-determined ode, 
consisting of its first ord[F] many equations, plus a purely algebraic equation, given by 
the square-free polynomial q. 
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Let q = q\ • • • q r be the decomposition of q as a product of irreducible factors in the 
polynomial ring K[£/, v]. Since the variables U are algebraically independent modulo the 
ideal (q) C K[U, v], we have deg„ qi > for all integer i. For each factor qi let pi be the 
ideal [j\ , . . . , f old [F] , Qi] in K{U, v}. 

Proposition 9 The ideal [F] is a radical quasi-regular differential ideal in ~K{U, v} and 
its minimal primes are pi, • • ■ , pr- 

Proof. Let us define B(U) '■= Y\bj(U). From the particular form of the polynomi- 
als F we observe that the ring ¥L{U, v }/\F] is isomorphic to a subring of the localization 
of ~K[U, v]/(q) at the polynomial Bq' . Since q is assumed square-free, the ring ~K[U,v]/(q) 
has no non-zero nilpotent elements and the same property remains true for any localization 
of it. Therefore the ideal [F] is radical. Similarly, each ring ~K{U, v}/pi is isomorphic to a 
subring of the localization of K[J7, v ]/(qi) at Bq', which is a domain since q[ is irreducible 
in ~K[U, v]; thus the ideals pi are prime. 

From the previous argument, we observe that the canonical map K[f7] — > f}/[£] 
is injective and no polynomial in K[U] is a zero divisor of the ring WL{U, v}/[F]. Since \F] is 
a radical ideal, it has only finitely many minimal (hence associated) prime ideals (see |45t 
Chapter 1, §16]) and so, none of these minimal primes contains a non-zero polynomial 
in K[U). 

Fix a minimal prime p of [F]. Since q is an element of [F] C p, there is an irreducible 
factor qi of q lying in the ideal p. Moreover, exactly one of these irreducible factors belongs 
to p since they are pairwise coprime in lK([/)[t>] (otherwise, Bezout's Identity would imply 
the existence of a non-zero polynomial in p n K[E7], leading to a contradiction). We will 
show now that pj is included in p. 

Since the total successive derivatives of the polynomials /i , • • • , / rd[F] belong to p 

and [F] is a subset of p, it suffices to prove that, for all integer j, the total derivatives qi ^ 
belong to p (note that p is not a priori necessarily a differential ideal). This can be 
done by recursion in j. For j = there is nothing to prove. Otherwise, for j > 0, we 
have that q^ +1 ^ = J2\h\=j+i hl^hA Q^ hl ^ ■ ■ ■ Qr^ hr ^ is in \F] C p, which implies by induction 
hypothesis that the expression qi ^ +1 ^ Yli^ qi is in p. Since Ili^io Qi i s n °t i n P> we 
conclude that qi ^ +l ^ is in p. 

Again, from the special form of the polynomials F_, it is easy to see that the system 
is quasi-regular at each minimal prime differential ideal pi and then, [F] is quasi-regular 
(see Definition [[]) . ■ 

The previous proposition ensures that the hypotheses of |10|. Section 2] are fulfilled. Hence, 
all the considerations concerning the differentiation index, the order and the Hilbert- 
Kolchin regularity explained there can be applied to our new DAE system S. In partic- 
ular we can compute the differentiation index of S at each minimal prime p as in [10\ 
Section 3.1]: 

Proposition 10 Let p be a minimal prime differential ideal containing [F]. Then the 
DAE system S has p- differentiation index 1. 

Proof. Following Proposition [9J let qi in K[C7, v] be an irreducible factor of q such 
that pj = [/i, . . . , /ord[F]) Qi]- Then, the fraction field K(pi) of K{U, v} /pi can be identified 
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in a natural way with the fraction field of the domain K[£7, v]/(qi); in particular K.[U] may 
be regarded subring of K(pj). 

Let 3i be the square Jacobian matrix of the polynomials /i, . . . , / or d[F]+l with respect 

to the variables U and v. Then 3i is the diagonal matrix: 



Jl 



Id 







c 



where C is a diagonal matrix with the elements bj(U)q' in the diagonal. Since bj(U)q' is 
non zero in the domain K[C7, v]/(qi) for all index j, we deduce that 3i has rank ord[F] 
over the field K(pj). 

Consider now ^2 the Jacobian matrix of the 2(ord[F] + 1) many polynomials F_, F_ with 
respect to the 2(ord[F] + 1) many variables U, v, U^ 2 \v . We have: 



&2 



31 





* ... * 


3i 


* • • • q' 



Therefore the rank of Z2 over IK(pj) is equal to 2ord[F] + 1 (recall that the last row and 
column of 2\ are both 0). 

Then, the relation dim K ( p .) ker*3i = 1 + ord[F] — rankjq Pi )(3i) = 1 hold and thus, the 
equality dim^.) ker*^2 = 2(1 + ord[F]) — rankj^p.)^) = 1 hold. Thus, from [TQl Defini- 
tions 5 & 9 and Theorem 8] we conclude that the pj-differentiation index of the system E 
equals 1. ■ 

In order to illustrate the notions above, let us now consider the classical pendulum example. 

Example 2 — Let £ be the DAE system arising from a variational problem describing 
the motion of a pendulum of length 1, where g is the gravitational constant and the 
unknown A is a Lagrange multiplier: 



£ 



y 



— Xx 

{2) -Xy + g 



x 2 + y 2 — 1 



0, 
0. 
0. 



(3) 



Using the notation F := {x^ — Xx, y^ — Xy + g,x 2 + y 2 — 1} to denote the system defin- 
ing E, we consider in the sequel the differential ideal [F] in the differential ring M{x, y, A}. 
The differentiation index of E is a = 4, as shown in |10|. Example 2]. We have, 



[F]nR (1) = pr^FnR (1) 



{x 2 + y 2 



1, yx' 



yxx 



y, xx + yy, x 2 + y 2 



yg + A, A - 3yg) 



and U := {x, x} is a common transcendence basis of the fields of fractions of R^ 1 ) / ([F] n R^ 1 )) 
and R (2) /([^] n R (2) ). In addition, it is not difficult to see that {x, x, X^ , X^\ \&\ X^ } 
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is a transcendence basis of Prac(R( e+<T Vpr (<T+1) i ? ) = Frac(R( 6 Vpr (5) -F), that is, we can 
take W := {X^\ X^\ X^\ X^} as in Proposition El 

We specialize the variables W in a point W € Q 4 in the generators of pr&F so that 
the conditions in Theorem 1181 below hold (note that, due to the structure of the system F 
and its successive derivatives, any specialization leads to a reduced regular sequence). 

In order to compute a parametric geometric resolution of [F] | yy , we consider the linear 
form v := y that defines the primitive element v. The minimal polynomial of this linear 
form is 

q = v 2 + x 2 - 1. 

Now, for each variable z £ {x^\ 2 < i < 6; yV\ 1 < j < 6; X^ h \ < h < 2}, we have poly- 
nomials b z (U) in M[U] and a z (U,v) in R[C7, v] such that 



6 2 (C/)^(C/,u)z-a 2 (L/,t;) 



vanishes over V := V([-F]|w)- For instance, 



b !J 


= 1, 


H 


bx 


= 1-x 2 , 


ax 


h 


= 1, 


a x 


x (2) 


= 1-x 2 , 


a x (2) 


yW 


= 1, 


a y (2) 


A( 2 ) 


= 1, 





— 2xx, 

-2x 2 v + 2g(l-x 2 ) 2 , 
—6gxx, 

-2xx 2 v + 2gx(l - x 2 ) 2 , 
—2gx 2 v — 2x 2 , 
6g 2 x 2 v — 6gx 2 . 



ui 



2\2 



If we denote u\ := x and u 2 '■= x, the system E from Definition [3] is 

2v(l 

By inverting ^pjj{U, v) = 2v modulo q(U,v), we get the simplified semi-explicit system 

(1 



- 2u\ii2 2 v — 2gu\{\ 

V 2 + Ux 2 - 1 



0, 
0, 
0. 



- UX 2 )U2 



v 2 + u 2 



U2, 

((1 - ux 2 )gv - u 2 2 )ux, 
1. 



-Ex 2H 



In the following section we will show how such a DAE system allows us to recover 
information about the original system E. 



4 Recovering solutions of E from solutions of E 

In this section, we will show that almost any solution of the system E introduced in 
Definition [3] can be lifted to a solution of the input system E. Moreover, we will prove 
that almost any solution of E may be recovered from a solution of E; more precisely, 
our main result states that there is a dense Zariski open set O of the variety of initial 



17 



conditions such that for any point X in O there exists a unique solution of £ with initial 
condition X that can be obtained by lifting a solution of E. 

For the sake of simplicity we assume that the ground differential field K is a subfield of 
the field of complex numbers C and the solutions of the involved systems are solutions in 
the classical sense. The arguments we will use can be easily extended to any differential 
subfield K of the field of rational complex functions C(t) by considering t as a new unknown 
variable and adding the equation t = 1. 

In order to lift solutions ip of £ to solutions ip of S, we start by introducing a dense 
Zariski open subset of the hypersurface T~L that defines suitable initial conditions deter- 
mining those solutions of £ that we will be able to lift. 

Let xft\ < j < e- 1, 1 < i < n, be a variable that does not belong to the set U. 
Let Ass (pr^ cr+1 )i ? |w; ) be the set of the associated primes of the radical ideal pr^ CT+1 ^i ? |yv; 
and let R p, where p runs over Ass (pr( CT+1 )F| w ),bethe primary decomposition of pr( CT+1 ).F| w 
(see Proposition U] and Notation [5]) . Then, for each component p (which is a prime ideal 
with dimp = card U) there exists an irreducible polynomial pjip in ~K[U, xft'] that lies in p. 

Let pji be the least common multiple of the polynomials (Pjip) pe Ass(pr( CT + 1 )_F1w)- Note 
that pji is the product of the polynomials Pjip-, without repeated factors: 

Pji = Y\_Pjip f° r some subset A of Ass (pr^ 1 ^ F\w ). (4) 

pGA 

Therefore, we have that pji is in pr( a+1 ^F\ w • Moreover, observe that no p in Ass (pr^ a+1 ^F\ w ) 
contains j^fe ■ In fact, the following relation holds: 

9 PH — \^ 9 Pjip TT n , 

peA pVp,P'eA 

Now if p is in A then L] pjip/ is not in p. Since all the remaining terms of the sum are 

Xl pVp 

multiples of pjip they lie in p, and we conclude that -^ijj is not in p. The same argument 
runs identically if p is not in A. 

Consider now the rational (not necessarily polynomial) map <!> : H — > V associated to 
the parametrization of V from H. Observe that for any polynomial p in K[X^ e +^\W] there 
exists a non negative integer h (depending onp) such that B(U) h ^(U,v) h p(Q(U,v)) is a 
polynomial in v], where B(U) is the polynomial introduced in the proof of Proposition 

El " 

Notation 11 Let h be a positive integer such that 

P0 (U,v) := B(U) h U(U,v) h J] Sl)^)) 

is a polynomial in K[U, v, x^)] divisible by B^. Note that, since ^jjj is not 

in p for each primary component p of pM^" 1 F\y\> , the set 

Go ■= U n {p / 0} 
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is a Zariski open set which is dense in the hypersurface %. Observe that Qq is included in 
the definition domain of the rational map <3? . 

Consider the projection 7Ti : A( 1+cr ) n+ord [ F ] — > A^ 1+e ^ n to the coordinates (x,x, . . . ,x^). 
Since the relation [F] n R^ e ^ = p^ a+1 ^ F\w n R^ holds (see Proposition [6]) , we conclude 
that tti(V) = V([F] n R( e )). In the sequel, we denote by S = V([F] n R (e) ). In particular, 
if the point (Uo, vq) is in Qq, the point -k\(&(Uq,vq)) is in S. 

Now we will show that an analytic solution <p of the DAE system E such that the 
point (f(0) is in Qq can be lifted to a solution (p of the DAE system E: 

Theorem 12 Suppose that ip : (-e,e) -»• C 1+ord [ F l is an analytic solution of the DAE 
system E stic/i i/iai is in £?o for all t and let 

V>{t) = W*), <Pe(t)) := vri o $ (£(t)) 

where <pj = (<Pj,i, ■ ■ ■ , Vj.n) for j = 0, . . . , e. TTien y?o : ( — £j s) — > C n is a weZZ defined 
analytic function that is a solution of the input DAE system E. 

Proof. Since <p(t) is in the subset Qq of Dom(<i>), the map ipo is well defined and analytic. 

In order to prove that tpo is a solution of E, first note that, since the image of tt\ o $ o </j 
is included in the variety 5, every polynomial in [F] n R^ e ^ vanishes at (7Ti o $ o ip)(t) for 
all i. In particular, this holds for each equations in E. 

Then, it suffices to show that the coordinate functions of tp verify the following relations 

= fj+iM ( 5 ) 

for every integer < j < e — 1, 1 < i < n and for any t £ (— e,e). To do so, one can 
consider two cases. 

First suppose that i,j are such that the variable x^ is an element, say Uh, of the 
transcendence basis U of Frac(R( e-1 - ) /[-F] D R^ 6-1 )) chosen in subsection 13. 11 Hence <Pj t i(t) 
is equal to iph(t) and relation ([5]) agrees with the equation corresponding to Uh in E after 
the specialization (U, v) i— > <^(i). Therefore this relation is satisfied because <p is a solution 
of E. 

Suppose now that Xi^ is not an element of U and let pji be its associated minimal 
polynomial in [F] nK[U,Xi^} C [F]nR( e ^ as in @. By taking the total derivative 
of pji, we obtain the following relation: 

Therefore ^pj j ((/?(£)) = for any t 6 (— e, e); more precisely, if ^ stand for the coordi- 
nates of (p(t) corresponding to U and U respectively, we have p'ji(<pu , y>j,i, = 0. 
Since we have already proved relation © for every u^, it follows that 

^ p ' ' d. Q p m ' 
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On the other hand, the polynomial p,j is in [F] n and, therefore, pji(ip(t)) = 0. By 
differentiating this identity with respect to t we obtain: 

Y.y^^u{t)^ 3 At))^u h {t) + ^^{<pu{t),wM)^jAt) = o- (7) 

h h % 

Since we assume that ip(t) is in Qq, in particular we have: 

Now, relation (J3J) is an immediate consequence of identities © and (f7|). ■ 

We have already shown above how we can recover a solution of the original system E from 
a solution of the new system E. Now we will show that almost every solution of £ can be 
recovered from a solution of E. We will apply the results of uniqueness and existence of 
solutions contained in Appendix IB1 

Recall that [F] is the differential ideal of K{£/, v} defined by the system E (see Defini- 
tion 1 page ED. Let 5 be the variety V([F] n K[?7, u, f/, u]) in A 2ord [ F l +2 and tt : 5 -> % 
the projection (J7, v, U, v) >->■ ({7, 

Theorem 13 There exist dense Zariski open sets O C S and O C S such that, for 
every point (Xq, . . . ,X e ) in O, there exist e > 0, a point (Uq,vq,Uq,vq) in O and an 
analytic function ip : (— e, e) — > C 1+ord ^ that is a solution of the system E with initial 
conditions ip(0) = (Uq,vq) satisfying: 

• {<£,'£) ■ (-e,e) ->■ Q, 

• (7Ti o $>)(ip{t)) = ((f(t), . . . , ip( e \t)) : (— e, e) — > O, where tp is the unique analytic so- 
lution of the system E with initial conditions (f(0), . . . , Lp( e ~ 1 ^ (0)) = (Ao, . . . , Af e _i). 

Proof. Let Q be the dense Zariski open set of regular points of S where the projection to 
V([F] n R^" 1 )) is unramified and Qq be the dense Zariski open set of regular points of S_ 
where the projection tt is unramified. 

Let us denote by O the dense Zariski open subset Qq D {po ^ 0} of S_, where pq denotes 
the polynomial in K[£/, v] introduced in Notation [11] page [THJ Note that, since all the 
points in 5 satisfy that B(U)^(U,v) ^ 0, the morphism <I? : T~L — > S is an isomorphism 
between the set tt(O) (denoted by Q\ in the sequel) and its image. Let us denote by O 
the set Q n {tt\ o <&)(Qi). We have then the following situation: 

V D $(Qi) A Qi C 7/ 

7T1 4 t 7? 

«S D O c 5 

Let Af be the point (Ab, . . . , X e ) in the set C Q. By Theorem [19] page [26] there 
exist a real e > 0, an open neighborhood 0^ of Af, Ox cOc5, and a unique analytic 
solution (fx ■ (— e, e) — > C n of E such that the image {(fx, ■ ■ ■ , Px^) e) is in Ox and 
the relation (<p x (0), . . . , ^ (e_1) (0)) = (X Q , . . .,X e -i) holds. 
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Since the point X is in tt\ o (^(Q-l), there exists a point £ in J 4(°"- e ) n + ord L F ] such 
that (Xq, . . . ,X e ,£) is in <3>(Qi) C V. Then, there is a point (Uo,vo) in Qi such that 
<J>(ZYo,^o) = (^0) • • • and, since Q\ = tt(0), the relation (Uq,vq) = 7t(Uq, vq,Uq, vq) 

holds for some point V = (Uo,vo,Uo,vo) in O. Recalling that O is a subset of Qq, by 
Remark [2 there exist a real e > 0, an open neighborhood Q-p C O of P and an ana- 
lytic solution ip : (—£,£) — > C 1+ord [ i? ] of (£) with initial conditions 99(0) = (Uq,vq) such 
that the image (92, <£>)(— £, £) is in Q-p. Then, for every t G (—§_,§), the point ip(t) is in 
the subset 7r(0) of C?0j where Go is the dense Zariski open subset of T~L from Notation [TT1 
page [TBI 

Now, Theorem 1121 implies that the relation tt± o = ((p(t), . . . , (p^ (i)) holds 

for t G (— e, e), where 93 : (— e, e) — >• C n is a solution of E. Since -k\ o $(99(0)) = tti o $(ZYo) ^0) = 
-^1(^0, • • • , Af e , ^) = X, which lies in O, taking a smaller e if necessary, we get that for 
every t G (-e,e) the point (<p(t), . . . , </? (e) (*)) 

is in 7Ti o ^(QjJ n Q = O. Moreover, 
as (<p(0), . . . , (p^ e ~^ (0)) = (Afo, . . . , X e ^i), if we denote min{e;, e} by e, using the uniqueness 
statement of Theorem PT9l we conclude that (fx(t) = (p(t) for every t G (— e, e). ■ 

Conclusion 

In this paper we presented a new index reduction method for a class of implicit DAE 
systems which is based on a characterization of the differentiation index from an algebraic 
point of view. We proved that any of these systems is generically equivalent to a first order 
semiexplicit DAE system with differentiation index 1 and we described a probabilistic 
algorithm to compute the index and the new DAE system by using the Kronecker solver 
for polynomial equations. 

Our results rely on some a priori hypotheses on the considered differential system, for 
example the primality of the ideal [F], which seems natural in practice, and the primality 
of the prolongation ideals, which is in general quite difficult to test. Assuming an admis- 
sible initial condition for our system to be known, polynomial time numerical methods of 
resolution, such as the one described in [8] or power series computations ([I]), may be used 
in order to obtain a solution. Since it is always possible to test, by simple substitution in 
the input equations, if a solution of the new system is actually a solution of the original 
DAE, one can attempt to use our method even if all our requirements on the considered 
system are not guaranteed. 

A next step in a future work would be to generalize the method to positive differential 
dimension and to regular components of systems, without any extra technical hypothesis. 

Urbs RomcE non uno die condita fuit. 

A On the specialization of free variables in a regular se- 
quence 

This appendix deals with Bertini-type results from Commutative Algebra that justify 
the random evaluation of suitable free variables made in Proposition 0] page [9l We have 
decided to include them for the sake of completeness and the lack of adequate references. 
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Throughout the appendix, K denotes a field of characteristic and the ring of 

polynomials in the variables X := {x±, . . . , x n } with coefficients in IK. 

We start recalling a well-known result concerning the behavior of radical ideals under 
field extensions of K (see [331 §27, (27.2), Lemma 2] and H3 Volume 2, Chapter 10, §11]): 



Theorem 14 Let IK C K* be a field extension and let I be a radical equidimensional 
ideal in K[X]. Then, the subset I* := J(8)K* of K*[X] zs a/so a radical equidimen- 
sional ideal. Furthermore, if I is a prime ideal and Y is a subset of the variables X 
such that K[Y] K[X]/I, then K*[Y] ^ K*[X]/p for any prime ideal p in Ass/* (in 
other words, the transcendence of the variables modulo I is preserved modulo any primary 
component of I*). 

Notation 15 For an ideal I C K[X], a subset W C X of cardinality s and any W G ~K S , 
we denote by 7|w the ideal ofK[X \ W] obtained after the specialization W i-> W in the 
polynomials of I. 



Lemma 2 — Let I be a radical equidimensional ideal in ¥L[X] of dimension s. Let W be 
a subset of s variables in X such that the canonical morphism K[W] — > ¥L[X]/I is injective. 
Then, there exists a nonempty K-definable Zariski open set O in A s such that /|w is a 
radical 0-dimensional ideal in K[X \ W] for any point W in On 1C S . j 

Proof. Let I denote the ideal / ® IK(VF) in K( W) [Y] where Y stands for the set of variables 
X \ W. From the hypotheses and Theorem [T4"l this ideal is a radical 0-dimensional ideal 
in the polynomial ring IC(Vl / )[y]. 

Then, the Shape Lemma implies that there exist a square-free polynomial q(W, v) 
in K[W][f], where v is a new single indeterminate, polynomials p y in K[W][w] for each 
variable y (zY, and a K-linear form £ in the variables Y such that the equality of ideals 

I=Lw,£),(^(W,£)y- Py (W,£)^ ) (8) 

holds in K(W)[Y]. Moreover, we may assume that the polynomials q, ^y — p y have triv- 
ial content in K[W]. From identity ([8|), it follows that there exists a nonzero denomina- 
tor g in K[W] such that the relation IK[X] g = (q(W, £), (^(W, £)y - p y (W, £)) y&Y ) K[X] g 
holds. 

Finally, let O be the Zariski open subset of A s where the product ( g discr„ q ) (W, £) 



is nonzero. Clearly for any point W in O we have 

I\ w = (q(W, £), (^(W, £)y - Py (W, £)) y& 

that is a 0-dimensional radical ideal in K[Y"]. ■ 

The previous Lemma can be generalized as follows: 

Lemma 3 — Let I be a radical equidimensional ideal in K[X] and W be a subset of s 
variables in X such that the natural morphism K[W] — > ¥L[X]/I is injective (in particular 



22 



the relation s < dim/ holds). Then, there exists a nonempty K-definable Zariski open 
set in A s such that, for any point W in D K s , the ideal J|w has dimension dim / — s 
and its primary components of maximal dimension are prime ideals. j 

Proof. Let us denote dim/ by d. If d = s, the result follows from Lemma [2j Sup- 
pose now that s < d. Let us denote by A the set of primary components p of / such 
that K[W] —> K[X]/p is a monomorphism. Since the canonical map K[W] —¥ M.[X]/I is 
assumed to be injective, the set A is a nonempty subset in Ass /. If we consider the 
canonical projection ir : A n — > A s the prime ideals in A correspond exactly with the com- 
ponents V(p) of V(I) such that ir(V(p)) = A s . Therefore, the set \Jp^ A ^(V(p)) * s a P ro P er 
closed subset of A s ; let Oq denote its complement in A s (that is a nonempty Zariski open 
set). 

Let U denote a set of d — s K-linear combinations of the variables X \ W such that the 
natural morphism K[J7, W] —> M[X]/p is a monomorphism for all p in A. After a change 
of variables we may suppose without loss of generality that U is a subset of X \ W. Let 
7 be the set of the n — d many remaining variables X \{U,W}. 

By Theorem ll4l the ideal / := / ® K([/) is a radical equidimensional ideal in M.(U) [W, Y] . 
Moreover, if p is in A, we have that p ® K(J7) is a prime ideal of dimension d — card U = s 
and so, / is also s-dimensional. Therefore, the ideal / C K(C)[W, Y] and the variables W 
meet the hypotheses of Lemma [2] and then, there exists a Zariski nonempty open set Q 
in A f. , definable over K(C7), such that /|w is a radical O-dimensional ideal in K(?7)[y] for 

every W in Q n K(C/) S . Let Oi be a Zariski nonempty open set in A s definable over IK such 
that 0\ is a subset of Q. Then, if W is a point in 0\ nl s , the primary components of /| w 
in M>[U, Y] with no nonzero polynomial pure in the variables U are (d — s)-dimensional 
prime ideals. 

On the other hand, since the projection ir : A n — > A s induces a dominant regular mor- 
phism of varieties tt : V (/)—>• A s , by the theorem on the dimension of the fibers (see for 
instance [50j Chapter I, §6.3, Theorem 7], suitably adapted to the case of a ground field 
not necessarily algebraically closed), there exists a K-definable Zariski nonempty open 
subset 02 of A s such that for every point W in 02 the fiber 7r _1 (W) is K-definable and 
geometrically equidimensional of dimension equal to dimV(I) — dimA s = d — s. 

Summarizing, if the point W is in O n 0\ n 2 n K s , then the ideal I + (W -W) 
satisfies: 

1. Any of its primary components of maximal dimension contains an isolated prime p 
in A (that is a subset of Ass /) such that K[U, W] ^ K[X]/p (since W is in O n 02 
and due to the choice of the variables U). 

2. Its primary components of maximal dimension which contain no nonzero polynomials 
in M.[U] are prime (because W is in 0\). 

3. It is geometrically equidimensional of dimension d — s (since W is in 02)- 

From condition (2), the Lemma will be proved if we are able to exhibit a "K-definable 
Zariski nonempty open set contained in Oq n 0\ D 02 such that, for any point W lying in 
this open set, all the isolated components of I + (W — W) contain no nonzero polynomials 
ofK[U]. The remaining part of the proof is devoted to showing this fact. 
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Let p be a prime ideal in A. We have the following injections K[W] ^ K[U, W] ^ IK[X]/p, 
and each variable yi in Y, with i = 1, . . . ,n — d, verifies a polynomial equation p%{yi) = 
modulo p, where pi is a nonzero element in K[J7, W][T] with deg T pi > (here T denotes 
a new variable). We denote by lc(pj) the leading coefficient of p, in ~K[U, W]. 

Claim: There exists a ^-definable nonempty Zariski open subset Oi$ in Oq n 0\ n 02 
suc/i £/ia£ /or any pomi W in 0, jP n and /or any irreducible component C o/7r _1 (W) n V(p) 
(necessarily of dimension d — s), we have that lc(pj) does not vanish identically over C. 

Let us prove this claim. To do so, consider the set V(p) n V(lc(pj)) in A n . If this set 
is empty there is no component of 7r _1 (W) n V(p) contained in V(\c{pi)) and the claim 
follows. On the other hand, if the set V(p) n V(\c{pi)) is nonempty, since we have the 
injection K[?7, W) ~K[X]/p, then lc(pj) is not a zero divisor modulo p and so this algebraic 
set must be an equidimensional algebraic variety of dimension diml/(p) — 1 = d — 1. We 
consider two cases: 

• If the relation 7r(V(p) n V{\c{pi))) = A s holds: from the theorem on the dimension 
of fibers applied to the restriction of the projection ir to V(p) n V (lc(pi)) , there ex- 
ists a K-definable nonempty Zariski open Oj jP in A s , that can be assumed contained 
in Oq Pi 0\ n 02, such that, for any point W in O^p C\K S , any component of the 
fiber 7r _1 (W) n V(p) n V(lc(pi)) has dimension dim V{p) n V(\c(pi)) - s = d- 1 - s. 

Now, if C is an irreducible component of 7r~ 1 (W) n V(p), since W is in O2 the 
relation dimC = d — s holds. If lc(pj) vanishes identically over C, then C is contained 
in 7r~ 1 (W) n V(p) n V(lc(pi)) that is a (d — 1 — s)-dimensional variety. So C is not 
contained in V(lc(pi)). 

• If 7r(V(p) n V{\c{pi))) is a proper subset of A s : define the set Oj jP to be the open 

Zariski set O n d D C 2 n vr(F(p) n y(lc( K ))) C . Clearly 7r -1 (W) n F(p) n V(lc(pi)) 
is empty and so Oj jP works. 

Hence, our Claim is proved. In order to finish the proof of the lemma, consider the 
nonempty Zariski open set O := f]i P eA ®i,P- ^ suffices to prove that for any W in O and 
any isolated primary component q of I + (W — W) C K[X], the relation K[U] n q = {0} 
holds. The primary ideal q defines an irreducible component C of the fiber 7r~ 1 (W) D V(I); 
therefore, C is a subset of V(p) for some p in A (since W is in Oq) and the relation dim q = d — s 
holds because W is in 02. Hence we have the sequence of natural morphisms 

K[U, W] ^ K[X]/p -> X[X]/Vq = K[C], 

where the last morphism is the projection to the quotient; in particular if we call (p the 
composition of the morphisms we have (j>(W) = W and the coordinate ring K[C] is gen- 
erated as a K-algebra by (j>(U) and the class of the variables Y. From the definition of O 
and the previous Claim, it follows that the class of each is algebraic over the sub- 
ring K[(^>({7)] of K[C] and so, the relation d — s = dimC < card U = d — s holds. Hence <j> 
is a monomorphism and in particular q does not contain polynomials in K[{7] \ {0}. ■ 
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Corollary 16 Let I be a radical equidimensional ideal in K.[X], and let U and W be 
two disjoint subsets of variables such that the natural morphism K[{7, W] — > M,[X]/I is 
injective. Then, there exists a nonempty ¥L-definable Zariski open set O in A s , where s 
is the cardinality of W , verifying the previous lemma and such that, for every point W 
in O n K s , the morphism K[U] ->■ K[X]/(I + (W — W)) is injective. 

Proof. From Theorem Q3] we are able to apply Lemma [3] to the field IK := K(£7) and the 
ideal I := 7 <8> K C ~K[X \ U\. Therefore we obtain a nonempty Zariski open subset O 
of A~. On the other hand, we may also apply Lemma [3] to the ideal 7 and the variables W 
over the ground field K, obtaining an open set Oq. 

Now take O an arbitrary K-definable nonempty Zariski open set contained in O n Oq. 
It suffices to see that K[U] is included in the ring K[X]/(I + (W - W)) for W in O, which 
is immediate from the fact that W is in O and, in particular, there exists a component 
of maximal dimension of 7 + (W — W) which contains no nonzero polynomial pure in the 
variables U. ■ 

If the ideal 7 is prime, the proof of Lemma [3] allows us to show a more precise version of 
the previous Corollary: 

Corollary 17 Let 7 be a prime ideal in K[X] and let U and W be two disjoint subsets of 
variables such that {U, W} is a transcendence basis of the fraction field o/K[X]/7. Then, 
there exists a nonempty ^-definable Zariski open set O C A s such that, for every point W 
in O n K s , the ideal 7|w has dimension dim 7 — s and its primary components of maximal 
dimension are prime ideals containing no nonzero polynomial pure in the variables U . 

Proof. Simply observe that in the proof of Lemma [3] the subset of associated primes A 
is the unitary set {7} if 7 is assumed to be a prime ideal. Therefore the Claim and the 
remaining part of the proof run mutatis mutandis. ■ 

Now, we can prove the main result of this appendix: 

Theorem 18 Let n > 2 and X a set ofn indeterminates over a field K. of characteristic 0. 
Let {/i,...,/ r } be a reduced regular sequence in M.[X] (that is, a regular sequence such 
that the ideals (/i, . . . , fi) in ~K[X] are radical for i = 1, . . . , r). Let W be a subset of s 
many variables in X, with s < n, such that the canonical map IK[W] — > . . . , f r ) 

is injective. Let us denote by Y the set of remaining variables X \ W . 

Then, there exists a nonempty ^K-definable Zariski open set O C A s such that, for 
every point W in O n K s and all i = 1, . . . ,r, the polynomials /i(W, Y), . . . , /i(W, Y) 
form a reduced regular sequence in the polynomial ring K[Y]. 

Proof. We prove this theorem by recurrence in r. 

If r = 1, since the polynomial f\ is assumed to be square- free, we take O as the 
projection of {d\scx y {fi) ^ 0} to A s , where y is any variable in Y which appears in f\. 

Assume the result holds for an integer r > 1. Let /i, . . . , f r +i be a regular sequence 
in ~K[X] such that the ideals (/i, . . . , fi) are radical for i = 1, . . . , r + 1 and let W be a 
subset of X such that the canonical morphism K[W] ^->- (/i, . . . , f r +i) is injective. 

In particular, K[W] > K[X]/(/i, . . . , f r ) is injective too. Hence, by the induction hypoth- 
esis, there exists a nonempty K-definable Zariski open subset 0\ of A s such that, for any 
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point V\ in 0\ and for any i = 1, . . . , r, the polynomials fi(Vi, Y), . . . , fi(V\,Y) form a 
regular sequence which generates a radical ideal in K[Y]. 

From Macaulay's unmixedness theorem (see for instance |27| Chapter VI, §3, Theo- 
rem 3.14]), the ideal I := (f\, . . . , f r +i) is equidimensional of dimension n — (r + 1) and 
so, the hypotheses of Lemma [3] are met for the ideal / and the variables W. Then, there 
exists a US-definable nonempty Zariski open subset 02 of A s , such that for any point V2 
in O2 n K s , the primary components of maximal dimension of (fi(T J 2, Y), ■ ■ ■ , fr+li^jY)) C 
are prime ideals of dimension dim! — s = n — (r + 1) — s. 

Take := 0\ n 02- We will show that this open set verifies the statement of the 
theorem for r + 1. Let W be a point in On K s . Since fi(W, Y), . . . , / r (W, Y) is a regular 
sequence that generates an equidimensional radical ideal (because W is in 0i), in order to 
prove that / r+ i(W, F) is not a zero divisor modulo (/i(W, Y), . . . , / r (W, Y)), it suffices to 
show that the dimension drops by 1 when this polynomial is added. This follows directly 
from the fact that W also belongs to 02- Therefore /i(W, Y), . . . , f r+ i(W, Y) is a regular 
sequence in K[Y]. In particular, the generated ideal is unmixed by Macaulay's theorem, 
and then it is radical since W E 02- ■ 



B Existence and uniqueness of solutions 



Several previous articles consider the problem of the existence and uniqueness of solutions 
for first order implicit DAE systems (see for instance [431 |4T)| |4"T1 110j). By adding new 
variables for the higher order derivatives in the usual way, these results can be extended 
to DAE systems of arbitrary order. For instance, let E be the DAE system introduced in 
Section 12.21 assuming that the hypotheses of Sections 12.11 and 12.41 are fulfilled and that K 
is a subfield of C. Then we have the following generalization of |10|. Theorem 24]: 

Theorem 19 Let Vq C A ne and V\ C A™( e+1 ) be the algebraic varieties defined by the 
ideals [F] n R( e_1 ) and [F] n R^ e ^ respectively, and let ir : Vi — > Vq be the projection to the 
first ne coordinates. 

Then, for every regular point X := (Xq, . . . , X e -i, X e ) in Vi where the projection ir 
is unramified, there exist e > 0, a relative open neighborhood C Vi of X and a unique 
analytic function if : (—e,e) — > C n which is a solution of £ with initial condition 

^(O),...,^- 1 )^)) = (Xo,...,Xe-l) 

such that {^{t),...,^ e ~ l \t),^ e \t)) is in for all t. 

Proof. We make a straightforward change of variables in order to obtain an equiv- 
alent first-order system: for each i G {0, ...,e — 1}, consider a new set Zi of vari 
ables (2^1, 
system 



,Zi tn ) representing the derivatives X®, and let T be the first-order DAE 



Zi-l, 



r := 



fi(Z,Z) 



0, 
0, 



1. 



,e-l, j 



1, 



(9) 



fn(Z, Z) 



0, 
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where Z denotes Zq, . . . , Z e _\. We apply now to this system the existence and uniqueness 
result in |1CH Theorem 24], which holds in the case e = 1. In order to do so, let us verify 
that the required assumptions hold. 

Denote by 21 the differential ideal associated with the system T and consider the 
map T : K{X} -> K{Z} defined by 

= (,-e +1) VT 

Note that T is an injection that maps [F] = [fi,---,f n ] to 21. For each differential 
polynomial / in the expression T(/) — (T(/))' belongs to the differential ideal 

[zij - Zilij ; 1 < i < e - 1, l<j<n). 

This implies that the relation 21 = T([i ? ]) + [zij — Zil\j ; 1 < i < e — 1,1 < j < n] holds. 
In particular, we have the isomorphism K{X}/[.F] ~ IK{Z}/21 and then, if [F] is a prime 
ideal, then so is 21. 

Moreover, we have the identities: 

ani[z] = T ([F] n R( e -^) , 

21 n K[Z, Z] = (z id - Zilxj ; 1 < i < e - 1, 1 < j < n) + T ([F] n R (e) ) . 

Therefore, if the polynomials F e _i in R^ e_1 ^ and F e in R^ e ^ are generators of [F] n R^ 6 " 1 ) 
and [F] n R^ e ^ respectively, then T(F e _i) and 

G := {^1,1 — ^0,1; ■ ■ ■ j 2l,n — 2o, n , ■ ■ ■ , 2 e -2,l ~ • • • , Ze-2,n ~ Z e -l,n} U T(F e ) 

are generators of 21 n K[Z] and 21 n K.[Z, Z\ respectively. 

In particular, the Jacobian submatrix f§ has the block form 

( -Id (e _i) n \ 

V T( J D x(e) F e ) j' 

Thus, if V(2l n K[Z, Z]) in A 2ne and F(2t D K[Z]) in A ne are the varieties defined by the 
specified ideals, and X is the point in V^(2l n K[Z, Z]) corresponding to X £ Vi, the 
unramifiedness of the projection 7r : V\ — > Vo at X is equivalent to the unramifiedness of 
the projection vr : 1/(21 n Z]\ ->■ F(2lnK[Z]) at ^. Similarly, the fact that Af is a 
regular point of Vi implies that is a regular point of V(2l n Z]). ■ 

Remark 2 — In the case of first order DAE systems, the existence and uniqueness of 
solutions as stated in [10t Theorem 24] can also be extended to the case when the ideal [F] 
is not prime, but the system £ is quasi-regular. The result follows as in the proof of (TOj 
Theorem 24] : if p denotes a minimal prime differential ideal of [F] , then p plays the same 
role as the ideal Q in that proof and we can take d = ord(p). j 
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